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The spatial distribution of electric current under magnetic field and the resultant orbital mag- 
netism have been studied for two-dimensional electrons under a harmonic confining potential 
V(r) = mus^r 1 /2 in various regimes of temperature and magnetic field, and the microscopic 
conditions for the validity of Landau diamagnetism are clarified. Under a weak magnetic field 
(lo c <ujq, lo c being a cyclotron frequency ) and at low temperature (T <Hujq), where the orbital 
magnetic moment fluctuates as a function of the field, the currents are irregularly distributed 
paramagnetically or diamagnetically inside the bulk region. As the temperature is raised under 
such a weak field, however, the currents in the bulk region are immediately reduced and finally 
there only remains the diamagnetic current flowing along the edge. At the same time, the usual 
Landau diamagnetism results for the total magnetic moment. The origin of this dramatic tem- 
perature dependence is seen to be in the multiple reflection of electron waves by the boundary 
confining potential, which becomes important once the coherence length of electrons gets longer 
than the system length. Under a stronger field ( lu c > wn ), on the other hand, the currents 
in the bulk region cause de Haas- van Alphcn effect at low temperature as T < fiuj c . As the 
temperature gets higher ( T > huj c ) under such a strong field, the bulk currents are reduced 
and the Landau diamagnetism by the edge current is recovered. 

KEYWORDS: orbital magnetism, confining potential, Landau diamagnetism, edge current, multiple reflection, 
mesoscopic system 



§1. Introduction 

Bohr-van Leeuwen's theorem tells us that the orbital magnetism does not appear in classical 
theory. SB The physical argument for this fact is that the diamagnetic current due to cyclotron 
orbits of electrons in the bulk region is perfectly cancelled by the paramagnetic current due to 
skipping orbits near the boundary. The orbital magnetism, which is then possible only in quantum 
mechanics, was originally derived based on the Landau levels.!) In the actual derivation, effects 
of the boundary have not been taken into account explicitly. These effects of boundary on the 

* E-mail: ishikawa@watson.phys.s.u-tokyo. ac.jp 



1 



orbital magnetism have been later investigated by many authors an Above 
all, Kubcft has applied the Wigner representation to an electron system under a magnetic field 




and shown that, if a confining potential is slowly varying in space compared to the electron wave 
length, the Landau diamagnetism results. However, this treatment is not valid at low temperature 
as stressed by Kubo and clear from the expansion parameters. The magnetic moment at T = , 



actually indicated that under a weak field and at low temperature ( /cgT < energy spacing by 
a confinement), the magnetic moment of the whole system shows a large fluctuation, and as the 
temperature is raised under such a weak field, fluctuations disappear and Landau diamagnetism 
is recovered. Hajdu and Shapiro© pointed out by studying the case of a groove with a width 
L that the temperature below which such fluctuations appear is such that fc^T ~ h/rt r where 
Ttr = L/vf, the time of flight for electrons at the Fermi energy across the system. Stimulated 
by these indications, we will study in this paper the spatial distribution of current in a confined 
system in order to understand more details of such a variety of orbital magnetism. From the 
investigations so far, the shape of confining potential, whether it is harmonic or hard wall, is 
expected not to affect qualitative aspects of the characteristic features of the orbital magnetism. 
Therefore we assume a harmonic potential which makes the analytical studies possible, and then 
clarify the relationship between the orbital magnetism and the spatial distribution of current in 
various regimes of temperature and magnetic field. 

The organization of this paper is as follows. In §|2], we introduce the model and summarize the 
general feature of the orbital magnetic moment. In the spatial distribution of current will be 
studied in various regions of temperature and magnetic field, and the microscopic conditions for 
the validity of Landau diamagnetism are clarified. Summary is given in §||. 

We take units ks = 1 in the following. 

§2. Magnetic Moment of Two-Dimensional Confined System 

In this section, we introduce our model and explain the general properties of the magnetic moment 
of the system. We consider two-dimensional electrons confined by an isotropic harmonic potential 
and the magnetic field is applied perpendicular to the system. For simplicity, we consider spinless 
electrons with the total number Nq and neglect Coulomb interactions between them. We assume 
No is large enough that the difference between a grand canonical ensemble and a canonical one 
is of no importance and rely on a grand canonical one to derive a magnetic moment. The word 
"magnetic moment" in this paper implies a magnetic moment of the whole system. 

The Hamiltonian is written as 





(2.1) 
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where p and r are two-dimensional vector, m is the electron mass and (— e) is the electron charge. 
The radius of the system, R, is classically defined as 



(2.2) 



where /i is the chemical potential. 

By taking a symmetric gauge A = x r, we can obtain an eigenfunction diagonal with respect 
to the angular momentum a as 
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where polar coordinates (r, 0) are used, and n = 0, 1, 2, ... , a = 0, ±1, ±2, ... and I = y/h/muj , to - 
a/o^T-P (2wo) 2 j w c = eH/mc being the cyclotron frequency and Ln is the Laguerre polynomial 
The eigenenergy of this state, E na , is given by 
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Especially under a extremely strong field (u> c 3> wn)> this eigenenergy E na becomes Tioj c (ji -\- 1/2) 
for negative a. So, the quantum number n corresponds to the Landau level index. 



By use of eq. (2^4), the thermodynamic potential Q is written as 
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where (5 = 1/T and \x is the chemical potential adjusted to fix the average electron number to iVo 
at each values of H and T. 

From the thermodynamic potential, the magnetic moment M is given as 
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where f(E) is the Fermi distribution function. 

Applying the Poisson summation formula© to the sum over n and a in eq. ( |2.5| ) (for details, see 
Appendix |~A| ) , one obtains 
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where w± 
in f2 



( fc + ^0 sinh(27r 2 ///3^ (T ) 
(a; ± co c )/2 and the relation S> T is assumed. In the derivation of oscillatory terms 



we used the approximation as follows, 
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where 6[x] is the Heaviside function. 

As is known from eq. (2.8), 0,q is dependent on a magnetic field only through the chemical 
potential n(H, T), which is regarded as a function of the field and temperature and almost a 
constant as a function of the field when fi/huj c 3> 1. Therefore the contribution of fio to the 
magnetic moment is much smaller than those due to Q,l and £l osc - 

Based on this result, the field dependence of the magnetic moment is classified into three re- 
gions; " Mesoscopic Fluctuation (MF)"," Landau Diamagnetism (LD) " and " de Haas-van Alphen 
(dHvA)", as is shown in Fig. 1. "MF" corresponds to the region as T < Tilu- , which implies 
T/hujQ <^ 1 under a weak field ( uj c /loo < 1 ) and T/Tiujq <^ (uj c /ljq)~ 1 under a strong field 
( lo c /loq > 1 ), while "LD" corresponds to the region as T > huj + , which requires T/Hloq > 1 under 
a weak field (u; c /ujq < 1 ) and T/Tiujq > uj c /ujq under a strong field (lo c /ujq > 1 ). The other region 
in Fig. 1 is "dHvA" . 

Fig. 2 shows field dependences of magnetic moment at various temperatures; (a),(b) and (c),(d) 
are respectively under a weak field ( u c < ujq) and a strong field ( cj c > ujq). Nq is set at 5000 
in this calculation. At first, we focus on the region of a weak field ( to c < ujq) shown in (a) and 
(b), where (a) corresponds to "MF" and (b) ranges from "MF" to "LD". At low temperature as 
T < Tiujq corresponding to "MF", the magnetic moment shows a large fluctuation with respect 
to the field, as Yoshioka and Fukuyama have shown. In this "MF" region, all the oscillatory 
terms in £l osc contribute to the magnetic moment leading to such a large fluctuation, in addition 
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to Oi = —1/2 • XlH 2 ( XL = — l/3 • Dq[j? b , the Landau diamagnetic susceptibility, where D = 
/i / (hu!o) 2 , the density of states at Fermi energy ), which leads to the usual Landau diamagnetism. 
Particularly at much lower temperature ( T <C Tiloq ) and under a much weaker field {oj c -C Uq ) , the 
magnetic moment shows a strong paramagnetism, which was first noticed by Meier and Wyder& and 
discussed by Budzin et al.& This strong paramagnetism is attributed to the rotational symmetry of 
the system. Hence, in the presence of weak disorder, a large spatial variation of magnetic moment, 
either paramagnetic or diamagnetic, is expected, which fact is the cause of a large variance of 
orbital susceptibility in the limit of weak magnetic field at low temperature. 0000® As the 
temperature is raised ( T > %ujq ) under such a weak field ( u c < loq), the fluctuation of magnetic 
moment is reduced and the magnetic moment becomes linearly dependent on the field, the slope 
of which gives the Landau diamagnetic susceptibilities xl corresponding to "LD" . In this "LD" 
region, the contribution of f2 osc are reduced and £Il becomes dominant. 

Fig. 2(c) and (d) show the field dependence of magnetic moment under a strong field ( lu c > loq) 
at various temperatures, (c) and (d) range from "MF" to "dHvA" and from "dHvA" to "LD", 
respectively. At low temperature as T < hto^ ~ Hluq/uJc, magnetic moment shows slow oscillation 
with a large amplitude and rapid oscillation with a small amplitude, as is shown in Fig. 2(c). The 
former oscillation with respect to the field is characterized by fi/huj ~ fi/huj c , which is caused by 
the periodic intersections of chemical potential \i by the Landau level ( the states characterized 
by n in eq. (|2.4|) ) and corresponds to the de Haas-van Alphen oscillation. The latter oscillation 
with respect to the field is governed by 4>/4>o where <fi is the total flux penetrating the system and 
(fro = hc/e, flux quantum. When the total magnetic flux is increased by (f>o, the degeneracy of each 
Landau level under Fermi energy is increased by unity. This causes the oscillation of the magnetic 
moment with a period </>o as a function of cj) at such a low temperature as T < huj- ~ Tiojq/oj c which 
is energy spacing between different angular momentum states at each Landau level. Physically, this 
oscillation is caused by a coherent motion of electrons along the edge and can be considered as the 
Aharonov-Bohm oscillation. Such a behavior under a strong field ( uj c > luo) has been noted at 
T = by Meir, Wohlman, and GefenEil As the temperature is raised so as huj- < T < huj + ~ ?luj c 
under such a strong field ( uj c > loq ) corresponding to "dHvA", the AB oscillation disappears 
and there only remains the dHvA oscillation. In this "dHvA" region, the first term in f2 osc has 
an appreciable contribution in addition to Ql- At much higher temperature ( T > huj + ~ hui c ), 
the dHvA oscillation disappears and magnetic moment shows a linear field dependence, i.e. the 
Landau diamagnetism as seen in Fig. 2(d) e.g. T /Tiloq = 4. 

§3. Spatial Distribution of Current and Magnetic Moment 

In this section, we study the relationship between a spatial distribution of current and a magnetic 
moment of the whole system, the latter of which is given by the thermodynamic potential in the 
previous section. 
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The spatial distribution of current in the system is given as follows, 

J(r) =Re/^t( r )Izf) fp( r ) + ^A(r)^) Mr) 

\ m \ c J 

= J e (r)e e , (3.1) 

where < • • • > denotes the thermal average and ip(r) is the field operator and eg = ^/\^\ as 
shown in Fig. 3. By use of the eigenstates in eq. (2.3) as the basis, Jg(r) is given as 



r \ 1 uj r / r 



a - + 
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where £ = ^h/muj , the characteristic length, and vq = ujq!; , the characteristic velocity of 
electrons. 

This local current J(r) induces a magnetic moment M{r) given as follows, 

M(r) = —r x J(r) 
2c 

2c 

= M z {r)e z . (3.3) 

Therefore, by integrating this magnetic moment with respect to r, we get the relation between the 
total magnetic moment M (z-component) and the local current Je{r) as 

M = J dS M z (r) 

= - dr r 2 J g (r). (3.4) 
c Jo 

This expression of magnetic moment derived from the local current density coincides with the one 



in eq. (2.6) derived from the thermodynamic potential (see Appendix B). 



3. 1 Region of Weak Magnetic Field ( oo c < luq ) 

Here, we focus on the properties of the spatial distribution of current in the weak field region 
( uj c <^ )• At low temperature ( T <^ hujQ, i.e. in "MF"), the magnetic moment shows a 
large fluctuation as a function of the field. The spatial distribution of current in such a situation is 
shown in Fig. 4(a), which indicates that Jg(r) can be either positive or negative, i.e. paramagnetic or 
diamagnetic, respectively. In the bulk region ( mainly, r < R), large currents flow paramagnetically 
or diamagnetically depending on r. This large bulk currents are sensitive to the strength of the field 
and lead to the large fluctuation of magnetic moment. This behavior of Je(r) is a characteristic 
feature of local currents in the region of "MF" . 

As the temperature is raised ( T > Tiujq ) under such a weak field, the fluctuation of magnetic 
moment disappears and the magnetic moment shows the Landau diamagnetism as discussed in 
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the previous section. According with this change of magnetic moment, the spatial distribution of 
current Jg{r) is changed as shown in Fig. 4(b); the fluctuating large bulk currents are immediately 
reduced and finally the diamagnetic current flowing along the edge (r ~ R) only survives. It is seen 
that this property of the current distribution in "LD" is characteristic of two-dimensional confined 
system as deduced from the Kubo's formula of the current distribution, which will be explained in 
the following. 

Based on the Wigner representation, Kubo& derived the analytic form of a current distribution 
proportional to a magnetic field and leading to the Landau diamagnetism in the system under a 
confining potential V(r), which is assumed to be slowly varying in space compared to the electron 
wavelength. The expansion parameters in this theory are ft 2 eH d ^ r - > /m 3 ^ 2 cT 5 / 2 , {fteH/mcT) 2 , 
ft 2 ( d ) /mT 3 and ft? d ^ /mT z . In the present model of harmonic confining potential, the 
parameters are explicitly given as oj c /u)q-{T/Tiujq)^/ 2 -R/^^ (uj c /ujq) 2 -(T '/ftu)o)~ 2 , (T/ftuio)~ 3 -(R/£) 2 
and (T/hoJo)~ 2 respectively, by replacing r with the system radius R. After all, T/ftujQ > (it!/£) 2 / 3 
and T/Tiojq > Uc/^o ar e required for the validity of the expansion in the Wigner representation, 
which are actually a part of the region "LD" in Fig. 1 ( However, such an expansion turns out 
to be an asymptotic one since our analytic result shows that ft, = is an essential singularity. 



See Appendix C for a detail. ). The spatial distribution of current, which is valid under such a 
condition, is given as follows in general, 

J(r) = ^cn 2 B Hx Vn(r), (3.5) 

where 

V(r)\f(E), (3.6) 



d is the dimension of the system and 7r corresponds to a physical momentum. The magnetic 
moment due to this current density becomes xlH, the Landau diamagnetism. 

In the present two-dimensional system, tv and E'-integration in eq. ( |3.6|) can be easily performed 
and n(r) is given by 

V/TTTfi 

n(r) = -^f(V(r)), (3.7) 
ft 

and J(r) becomes 

J{r) = \c&~-Hx W(r) • f'(V(r)). (3.8) 

Due to the factor f'(V(r)) which reflects Fermi degeneracy, this form of a current distribution 
means that the diamagnetic current flows only in the region V{r) ~ fi , i.e. the edge of the system, 
even if the potential is varying spatially in the bulk region. 

Applying the above formula to the present model V{r) = raio^r 2 / 2, we get J(r) as 
j, ^ 30 " C fry 1 r ePW-ri 
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where 
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(3.10) 



and jo = ewo/47r£ the characteristic current density. The current distribution calculated by eq. ( |3.9[) 
coincides with the one obtained from eq. (3.2). ( For example, at co c /ujq = 0.1 and T/Tiluq = 10, 
the agreement is within 0.5% numerically.) 

In comparison, the case of d = 3 which Kubo considered appears somewhat different, since the 
current both at the surface region ( the region as V(r) ~ [i ) and the bulk region can contribute 



to the magnetic moment. To see this, we approximate a Fermi distribution function in eq. (3.6) by 
the linearized form around the Fermi energy as 
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We assume the confining potential V(r) is only dependent on r in cylindrical coordinates (r,6,z). 
Then the current distribution is given by 
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where 7 = v / 27rw- 3 ^ 2 C/U^/3/i 3 and g(V(r)), the function corresponding to f'(V(r)) in the case of 
d = 2 is given as 
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In the bulk region where V(r) <C fi— 2T , g{V{r)) becomes 

2 



ff (V(r)) ~ -=%= « 1 , (3.14) 

while at the surface region, |V(r) — jx\ < 2T , g(V(r)) ~ 1 . Therefore, it can be said that 
the current causing the Landau diamagnetism is mainly induced at the surface also in a three- 
dimensional system, which is not so clear as the case of d = 2. 

At low temperature as T < Tllvq , on the other hand, the contributions of higher order terms 
in magnetic field become larger and the above formula of current density is not valid. Then, the 
current distribution changes dramatically at around the temperature T = hujQ under a weak field 
(uj c < ujq) as is shown in Fig. 4(a) and (b) and this change is clearly reflected in a magnetic moment 
of the system. Hajdu and ShapiroE!) studied a two-dimensional system under a confining potential 
V(r) = mu>QX 2 /2 ( i.e. harmonic groove ) with a width L x and an arbitrary long length L y by 
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imposing a periodic boundary condition in the y-direction, and pointed out that the temperature 
as T = hujQ below which magnetic moment shows a large fluctuation under a weak field ( uj c < uq ) 
corresponds to Ti/rtr where rj r = L x /vp, time of propagation for electrons at the Fermi energy 
across the groove. The physical reason why the field dependence of magnetic moment shows such 
a dramatic difference depending on the temperature is understood as follow. When the effect of a 
confining potential V(r) is considered perturbative, the thermal Green's function of a degenerate 
electrons has the following damping factor, 

G(x,t = 0) oc e~^ x . (3.15) 

Therefore, the length Uvf/kT = l c is regarded as the coherence length of degenerate electrons. 
Under the condition l c > L, where L is a system length, electrons near the Fermi surface can 
propagate from one side of the system to the other side with a small damping and experience the 
multiple reflection by the boundary potential as shown in Fig. 5(a). Therefore electrons near the 
Fermi surface, which play an essential role in the orbital magnetism of degenerate electron system, 
are strongly affected by the boundary potential. In this model, the condition l c > L corresponds 
to T < hujQ i.e. "MF". As shown in Fig. 4 (a), it is seen that this multiple reflection induces large 
currents irregularly distributed paramagnetically or diamagnetically in the bulk region, which are 
sensitive to the strength of a magnetic field and causes the large fluctuation of a magnetic moment 
as a function of the field. On the other hand, under the condition l c < L ( T > Twqq ), the effect 
of the multiple reflection by the boundary potential wall is reduced as shown in Fig. 5(b), and 
this is considered to lead to the suppression of the bulk currents and the recovery of the Landau 
diamagnetism. Considering the case of a harmonic groove system based on this idea, which Hajdu 
and Shapiro^ studied, it is natural that the relative magnitude of l c and the width L x affects 
the field dependence of magnetic moment while L y does not affect it qualitatively, since the re- 
direction is actually confined by a harmonic potential and a multiple reflection can happen only in 
the x-direction. 

Robnik§) discussed the size effect on the zero-field susceptibility in a system confined by a hard 
wall, based on the Green's function method, and concluded that at T = the contribution of 
the boundary wall to the susceptibility is always paramagnetic but with the order of magnitude 
(kpL)~ 1 compared to the Landau susceptibility, where kp is a Fermi wave number and L is a system 
length. As is clear from Fig. 2(a), this evaluation of a size effect due to the confining potential 
is too small. This disagreement can be attributed to the fact that the multiple reflection by the 
boundary potential wall is disregarded in ref. 9. 

3.2 Region of Strong Magnetic Field (lo c >ujq) 
Under a strong magnetic field ( uj c > ojq ), the energy spectrum in eq. fl2,4| ) becomes close to the 



one without the confining potential, i.e. the Landau level ( the states characterized by n in eq. (2.4) 
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). At low temperature as T < Tilo c , magnetic moment is affected by this Landau quantization and 
oscillates as a function of the field as is shown in Fig. 2(c) and (d) (This is the familiar de Haas-van 
Alphen effect). In this "dHvA", the spatial distribution of current is as is shown in Fig. 4(c). 
Several diamagnetic peaks can be seen and paramagnetic currents flow between them. These bulk 
currents are much larger than the edge current in "LD", although the diamagnetic currents are 
almost cancelled by the paramagnetic ones as an average. Each diamagnetic peak comes from the 



state with the same n in eq. (2.4) which corresponds to Landau level. In this model, each Landau 
level has a different degeneracy with respect to the angular momentum, a, due to the confining 
potential and the lower Landau level has the higher degeneracy. Therefore, the spatial extent of the 
lower Landau level with such a degeneracy becomes larger and is reflected in the current distribution 
as the spatially separated diamagnetic peaks. As the field is raised, the number of Landau level 
below the Fermi energy decreases and the diamagnetic peaks closest to the center of the system 
disappear one by one. This change of the current distribution causes the large oscillation of the 
magnetic moment as a function of the field. At much lower temperature in "MF" ( T < Tvuj-. ) 
where the AB effect appears, the current distribution is almost same as in "dHvA" region. 

At higher temperature as T > hiv c , the diamagnetic peaks in the bulk region are smeared out 
and there remain only edge current as is shown in Fig. 4(d). This edge current is spatially broad- 
ened compared to the one at lower temperature in "LD" shown in Fig. 4(b). This temperature 
dependence of edge current is understood from eq. ( |3.9[ ) where the spatial extent of edge current 
corresponds to | V(r) — fj, \ <T. 

§4. Summary 

We have studied the relation between the spatial distribution of current and the magnetic moment 
of the whole system in a two-dimensional electron system under an isotropic harmonic potential as 
V(r) = mwQf 2 /2. It is found that characteristic dependences of magnetic moment on temperature 
and magnetic field are clearly reflected in the spatial distribution of current. 

Under a weak field ( uj c < ujq), the field dependence of magnetic moment dramatically changes 
at around the temperature T = Hujq . In the low temperature region ( T <^ Twjq ), magnetic 
moment shows a large fluctuation as a function of the field, as was indicated by Yoshioka and 
Fukuyama.0 We attribute this large fluctuation of magnetic moment to a multiple reflection by 
the boundary confining potential, which becomes important once the coherence length of degenerate 
electrons l c = Tivf/t^T gets longer than a system length. In the present model, the system length 
is characterized by the radius R and the above condition leads to l c > R implying T < Tiujq . 
At such low temperature, it is seen that the multiple reflection induces large currents irregularly 
distributed paramagnetically or diamagnetically in the bulk region, which are sensitive to the 
strength of the field and cause a large fluctuation of the magnetic moment. As the temperature 
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is raised ( T > Tvujq ), the fluctuations of the magnetic moment are reduced and the usual Landau 
diamagnetism is recovered. Corresponding to this change, it is found that the large currents in the 
bulk region are immediately reduced and finally the diamagnetic current flowing along the edge ( 
in the region satisfying V(r) — // ) only survives, which leads to the Landau diamagnetism. This 
edge current is regarded characteristic of a two-dimensional confined system, as inferred from the 
Kubci's formula of the diamagnetic current distribution which is derived for high temperature as 
T > Tiujq ■ (i?/£) 2 ' 3 and hco c , where £ = y/h/muiQ , the characteristic length. From this formula, 
the cancellation of the bulk currents is seen to be due to the Fermi degeneracy. It is also noted 
that the persistence of the diamagnetic current at the surface is also seen in a three-dimensional 
system. 

Under a strong field ( oj c > ojq ), dHvA effect appears at low temperature as T < Tiuj c . In this 
situation, the current is distributed in the bulk region and diamagnetic peaks appear, although 
paramagnetic currents flow between the peaks, which almost cancel the contribution of the dia- 
magnetic peaks as an average. Here, each peak is due to the Landau level with different degeneracy 
because of the confining potential. As the field is increased, this diamagnetic peak disappears one 
by one corresponding to the decrease of the number of Landau levels under the Fermi energy, and 
this causes the large oscillation of magnetic moment (dHvA effect). At much lower temperature 
( T < %U— ~ Hluq/uJc , the energy spacing between different angular momentum states at each 



Landau level ), the small but rapid oscillation with a period </>o = hc/e appears as a function of 



at T = 0. This oscillation is caused by a coherent motion of electrons along the edge and can be 
called AB oscillation. On the other hand, at high temperature as T > hui c under such a strong 
field, the bulk currents are quickly reduced and the Landau diamagnetism by the edge current is 
recovered as in the case under a weak field. It should be noticed that this edge current does not 
cause AB oscillations, because the current is due to incoherent motions of electrons. 
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Appendix A: Derivation of Magnetic Moment from Thermodynamical Potential Q 

The thermodynamic potential f2 is given by 




the total flux cf) in addition to the dHvA oscillation, , which was ori; 



ginally found by Meir et a 
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where E na = hu{n + 1/2) + twj\a\/2 + fou> c a/2. 
By applying the Poisson summation formula 
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n+- = / dxF(x) + 2ReV(-l) fc / dxF( 
^0 V 2 ^ 7o *=! Jo 

to the sum over n, $7 is transformed as 



2irikx 



Q = - 



1 Z" 00 



Jo 

1 /"°° 



de 
de 



l + 2Re]T(-l) fc e 2wi/£ >fc 

fe=i 

00 

l + 2Re^(-l) fc e 2wifc 
fc=i 



I3Tluj Jo 
where to± = (lo ± lo c )/2. 

To the sum over a, we again use the Poisson summation formula 



log [ 1 + e -^ £ -^" 
00 

^ log 1 + e -^(^"+e~M) 
ct=± a=l 



as 



1 °° />00 00 />00 

-F(0) + V F(a) = / ds^s) + 2ReV / dxF(x)e 2rik . 

2 ^1 Jo ^ 



Then O is given as 

^ /■oo roo 



n = - 



roo roo 

- E ft Re E / de / dr? log 1 + e^^-") 



• log [ 1 + e - /3(e+,? -^ 



4 00 roo 00 /-oo 

— — Re Y(-l) k / dee 2 ^ ifc ^ -ReV / dry log 
flU}(j ^ Jo ^Jo 

We perform integrations over e and rj twice by parts, 



E 



1 _|_ e -y8(£+^-M) 



roo _ , 

/ de log 1 + e^ 4 -^ 
Jo L J 



■ e n ^ 



= _^i og [l + e -^-M) 
27ri/c L 



+ " Is* J /( *» 



flLO \ 

2irk 



e nui 



roo 



■/?&*-»?) (e* + l) 2 



■ 2irk e 
. e /3ft" ? . 



In the last integration term, we assume fj,^$> T and approximate as 

f°° 

\ ^ TTT T\2 • e f 



27T 2 fc 



sinh 



/ 27T 2 fc \ 
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where 0[x] is the Heaviside function. 

By using this approximation and the relations 



E 

k=l 



(-1) 



fc+1 



k 2 



TT y. J. /( 



1 



k=l 



k 2 



IT 

~6 



we get the form represented in eq. ( |2.7| ). 

Appendix B: Derivation of Magnetic Moment from Current Distribution 

By replacing r 2 /2l 2 in eq. (3.4) with p, magnetic moment is given by 



oo oo 



M 



■^E E 



n 



n=0 a=— oo 

p-integrations can be performed as 



( n + \a\ 



f(E n 



dp 



0J C 

a H p 



-Pp^L^pf 



' dp e"V |a| L^[P] Lt l) iP) = (n + H)! «- 



n ! 



and 



dp e- p pl a l +1 L(l a l)[p] 2 



dp e -v o|+i [^ a|+i) w-^n +1, w 

n + |a| + 1 ) ! ( n + |a| ) ! 



+ 



n l 



n + a 



n 



n ! 



- ■ (2n+\a\ + l) . 



Therefore, the magnetic moment becomes 

M : 



oo oo 



E 



wE E 

n=0 ct=— oo 
3E n 



a + — ■ (2n + \a\ + 1 



dH 



f{E na ) 



(B-l) 



(B-2) 



(B-3) 



(B-4) 



This coincides with the one derived from the thermodynamical potential in eq. (2.6). 

Appendix C: Wigner Representation for Landau Diamagnetism 

In the present model of the harmonic confining potential, the parameter region, where the 
Kubo's formula of current distribution leading to the Landau diamagnetism^ is valid, is T/hujQ > 
^ (R/t;) 2 ^ 3 = (2 • p/huJo) 1 / 3 and T/Tlluq > uj c /loq , which is dependent on the chemical potential, 
|U, i.e. the number of electrons. This region is only a small part of "LD" region ( T/Hujq > 1 and 
T/huo > oj c /ojq ) of Fig. 1. The reason of this apparent discrepancy is as follows. 

By differentiating the trigonometric functions in eq. (2.10) with respect to the field, it is found that 
VL osc contributes to the magnetic moment in the order of magnitude (2ir 2 T /Tiuj) / sinh(2-7r 2 T '/ 'huj) or 
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(27r 2 T/huj±)/ sinh(27r 2 T /hu±) relative to Oj,, and based on this fact the phase diagram in Fig. 1 is 
decided. These factors (2vr 2 T jfuS) / sinh(27r 2 T/ '%uj) and (2vr 2 T//ia; ± )/sinh(27r 2 T//ia; ± ) can not be 
expanded in terms ofhu/T and Hu±/T, which are seen to correspond to the expansion parameters 
in the treatment of Wigner representation. This implies that one cannot conclude whether the 
bulk magnetic moment shows the Landau diamagnetism or not, unless the contributions of higher 
order terms in h in Wigner representation are summed up to the infinite order. Actually, in the 
parameter region, 1 < T/huj < (i?/£) 2 / 3 and T/Hluq > lo c /loq, some expansion parameters in ref. 4 
get larger than 1, but even in such a case, the Landau diamagnetism is realized and diamagnetic 
current is mainly induced at the edge as seen in Fig. 4(b). 
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Fig. 1. Phase diagram showing characteristic regions of the field dependence of magnetic moment at various tem- 
peratures. 



Fig. 2. Magnetic field dependences of magnetic moment at various temperatures. The dependences under a weak 
field ( LJc < ioo ) are shown in (a) and (b), where (a) corresponds to "MF" and (b) ranges from "MF" to "LD". On 
the other hand, the dependences under a strong field (lo c >loo ) are shown in (c) and (d), where (c) and (d) range 
from "MF" to "dHvA" and from "dHvA" to "LD", respectively. The number of electrons is fixed to 5000. 



Fig. 3. A schematic illustration of the current flowing in the system. Je{r) is defined positive in the direction as 
here. 



Fig. 4. Spatial distribution of current at various temperatures. Here, (a) and (b) are the current distribution under 
a weak field (u) c /uio = 0.1), where (a) corresponds to "MF" and (b) ranges from "MF" to "LD". (c) and (d) are the 
current distribution under a strong field (ui c /u)o = 20 ), where (c) and (d) range from "MF" to "dHvA" and from 
"dHvA" to "LD", respectively. The number of electrons is fixed to 5000. Je(r) is normalized by jo = eu>o/47r£. 



Fig. 5. A schematic illustration of the reflection at the boundary confining potential; (a) l c (= %vf /kT) > L where 
the multiple reflection becomes dominant, and (b) l c < L where the boundary affects the behaviors of electrons 
only around the boundary. 
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